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ABSTRACT 

We present detailed chemical element abundance ratios of 17 elements with 8< Z <60 in three 
metal poor stars in the Ursa Minor dwarf spheroidal galaxy, which we combine with extant 
data from the literature to assess the predictions of a novel suite of galaxy chemical evolution 
models. The spectroscopic data were obtained with the Keck/HIRES instrument and revealed 
low metallicities of [Fe/H]=-2.12, -2.13 and -2.67 dex. While the most metal poor star 
in our sample shows an overabundance of [Mn/Fe] and other Fe-peak elements, our overall 
findings are in agreement with previous studies of this galaxy: elevated values of the [a/Fe] 
ratios that are similar to, or only slightly lower than, the halo values but with SN la enrichment 
at very low metallicity, as well as an enhancement of the ratio of first to second peak neutron 
capture elements [Y/Ba] with decreasing metallicity. The chemical evolution models which 
were tailored to reproduce the metallicity distribution function of the dSph, indicate that UMi 
had an extended star formation which lasted nearly 5 Gyr with low efficiency and are able to 
explain the [Y/Ba] enhancement at low metallicity for the first time. In particular, we show 
that the present day lack of gas is probably due to continuous loss of gas from the system, 
which we model as winds. 

Key words: galaxies: individual (UMi I dSph)-galaxies: dwarf-galaxies: chemistry and star 
formation history—Local Group—galaxies: abundances—galaxies: evolution 


1 INTRODUCTION 


The increasing number of detailed spectroscopic surveys inves- 


(dSph) galaxies (see e.^.. To 

stov et al. I 2 OO 4 I: iHelmi et al.l l2006l; 

iKoch et al.ll2006l: iMartin et al. 

I 2 OO 7 I: Walker et al.ll2009ab nrovide 


new opportunities to study their local formation environment in de¬ 
tail and the formation and evolution of galaxies in general. Despite 
their low luminosities, their gravitational potentials had to be large 
enough to confine their gas for sufficient time to allow an extended 
period of star formation (SF) to occur. This is indicated by the large 
spreads in their metallicities a nd by the presence of stellar pop¬ 
ulations as young as 1-2 Gyr (iTolstov et al. 2009|. and references 
therei n), even down to 200 Myr ( Stetson et al.l 19981 : Ide Boer et al.l 
I2OI2I) . Samples of ev en a few stars in faint dSphs can pro vide in¬ 
teresting information (iKoch et al.ll2008l : [Simon et al.ll201Cll) as they 


emphasise the diversity of the chemical inventories of individual 
dSp hs and guide us thro ugh the processes that affected their evolu¬ 
tion (iTolstov et al.ll2009l) . 


One of the most important observed signatures is the spread in the 
ratio of [cr/Fe] in dSphs where the downturn of [of/Fe] with in¬ 
creasing metallicity marks the time scale on which SN la started 
to contribute to the enrichment of the interstellar medium (I SM). 
Other mechanisms such as episodic SF (iHendricks et al.ll2014h or a 
change in the Initi al Mass Function (IMF) where the very massive 
stars do not occur (IMc W illiam et al .112013 ) might also be responsi¬ 
ble for the spread found in the abundances of the cr-elements in the 
present day dSphs. 


Although many isolated dwarf galaxies seem to be dom- 
inate d by exclusively old _ ^d/or young stellar popula¬ 

tions teemtez Llambav et al.ll2014h . among the dSphs, Ursa Minor 
(UMi) is a rare system which has exclusively old stars and hence a 
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Target 

I 

V 

K 

a J2000 

d (J2000) 

Teff (phot) 

Teff (spec) 

logg 

Vmic 

[Fe/H] 

Vr (kms D 

UMI718 

16.30 

17.46 

14.967 

15:02:25.58 

66:55:25.39 

4563 

4540 

1.1 

2.25 

-2.12 

-247.6 

UMI396 

15.59 

16.94 

14.018 

15:05:52.04 

66:41:30.05 

4312 

4290 

0.5 

2.35 

-2.13 

-228 

UMI446 

16.96 

18.07 

15.131 

15:02:01.53 

66:49:22.36 

4640 

4580 

1.6 

2.80 

-2.67 

-246 


Table 1. Observed stars: The K magnitudes are obtained from 2MASS photometry while I and V are from iKlevna et alJ jl998li . The preliminary estimates for 
the photometric temperature T(phot) from (V-I) and the surface gravity are calculated assuming a metallicity of [Fe/H]= -2 dex, mass of O.SMq and distance 
of 68 kpc.The T(spec), log g, Vmic and [Fe/H] are the values found by iteration of the best atmospheric models used to calculate the spectroscopic abundances 
(see text for details). 


candidate for being a dSp h similar to those wh ich contributed to the 
growth of the outer halo. ICarrera et ahl (l2002h suggested that 95% 
of its stars are older than 10 Gyrs and llkuta & Arimot^ (l2002l) 
constrain the duration of SF in UMi to be between 3.9 to 6.5 Gyr. 
In this paper, we add the chemical abundances of three red giant 
br anch stars in UMi to the eighteen stars from the previous studies 


bvICohen & Huang (2010l); Shetrone et al.l (l200ll):ISadakane et al.l 

(12004): Kirbv & Cohe 

rJ (I2OI2I). The chemical evolution models 

tested bv Kirbv et al.! ( 

201 ll) andiKirbv et al.! (I20131) showed that 


the metallicity distribution function (MDF) of UMi can only be 
explained with models without extended infall of gas. We present 
detailed chemical evolution models describing the nucleosynthesis 
of 12 elements for the combination of our new data set with the 
extant data, examine these findings and study the implications for 
individual elements. 

The paper is organised as follows: In Section |2] both the observa¬ 
tions and the methods used to obtain the chemical abundances are 
explained and the final abundances for or, Fe-peak and neutron cap¬ 
ture elements are given for three new stars. In Section[3l we explain 
the details of our chemical evolution models and compare them to 
the combined data set. In Section!?! we discuss our results and open 
questions. 


A E.P. EW [mA] 

[A] [eV] UMI396 UMI718 UMI446 


0 

6300.31 

0.00 

-9.780 

34.6 

- 

- 

Mg 

4571.10 

0.00 

-5.623 

- 

123.6 

93.5 

Mg 

4702.99 

4.34 

-0.440 

98.9 

111.2 

91.1 

Mg 

5528.41 

4.34 

-0.498 

132.4 

107.8 

- 

Mg 

5711.09 

4.34 

-1.724 

36.3 

- 

- 

Si 

7423.50 

5.62 

-0.314 

30.8 

- 

- 

Ca 

4435.69 

1.89 

-0.519 

- 

- 

72.7 

Ca 

4454.78 

1.90 

0.260 

- 

116.0 

- 

Table 2. 

Equivalent width measurements. Columns are: (I) Element; (2) 

Wavelength of the absorption line; (3) Excitation energy of the lower energy 


level; (4) Oscillator strength; (5,6,7) Equivalent width of the line measured 
for each star. This table is available in its entirety in the electronic version 
of the Journal. 


The continuum fit to each echelle order was obtained by dividing by 
high-order polynomials within the continuum task in IRAl(|. Sub¬ 
sequently, the velocity correction from vacuum to air velocities was 
made using the disptrans utility and the Doppler shift was com¬ 
puted from the average wavelength shift of a few of the strongest 
lines in several orders. These velocities are given in Table [T] 


2 OBSERVATIONS AND ANALYSIS 
2.1 Observations and data reduction 

The observations wer e made using the High Resolution Echelle 
Spectrograph (HIRES: IVogt et al.l(ll994h ) on the Keck telescope on 
the night of UT 2006-06-05. The t argets were selected fro m legacy 
KPNO 4 m MOSAIC imaging as in iwilkinson et al.l(l2004h . Details 
of the targets are given in Table [T] together with some of their basic 
characteristics. 

We used HIRES configured with the red collimator, the B2 decker 
set for a fully open slit, the kv408 filter with the cut off at 
450 nm, an echelle angle of OaU^.11197 and cross-disperser angle 
of OaU^.65899. We obtained a full wavelength coverage of 440- 
750 nm with a 10 nm wide gap at 590 nm. The exposure times were 
4000s for UMI718, 3600s for UMI396 and 3000s for UMI446. 

The raw data were reduced using the MAKEE data reduction pack- 
ag4E which performs standard steps such as bias and fiat field cor¬ 
rections, wavelength calibration via Th-Ar lamps, optimal subtrac¬ 
tion, and sky subtraction. This procedure resulted in signal-to-noise 
(S/N) ratios measured in the Ka order around 660 nm of 25, 30, 12 
per pixel, respectively, for UMI718, UMI396 and UMI446. 


2.2 Stellar parameters and abundance calculations 

We derived chemical element abundances via a standard equivalent 
widths (EWs) analysis t hat closely follow s the procedures used in 
our previous works (e.g. lKoch et al.|[200^ . All analyses e mployed 
the 2010 version of the stellar line analysis code MOOG (ISnedenI 
Il973[) . The EWs of the absorption lines were measured within the 
splot task in IRAE assuming a Gaussia n line profile. The line list 
was assembled from iMoore et al.l (1 19661) combined with the log gf 
values from the Kurucz database. The EWs for all the lines used 
in our analysis are given in Table |2j In order to avoid saturation or 
spurious, weak lines, we conservatively excluded all Ee I and Ee II 
lines with EW^ 160mA and EW<40mA from our analysis. 
Although we do not expect the effect of hyperfine splitting to be 
large in our spectra, they were investigated in the derivation of 
abundance ratios for Sc, Mn, Cu, and Ba, using wavelengths and 
log §f values taken from the Kuru cz database with the ex cep- 
tions of the Sc II line at 467.0 from iMcWilliam et"^ (1 1995!) and 
565.7, and 568.4 nm, calculated with the coefficients from jjohnsonl 
(l2002h . Except for the larger corrections for [Mn/Ee], which ranged 
from 0.09 to 0.15 dex, the inclusion of hyperfine splitting causes 


^ MAKEE was developed by T. A. Barlow specifically for 
reduction of Keck HIRES data. It is freely available on 
the World Wide Web at the Keck Observatory home page, 
http://www2.keck.hawaii.edu/inst/hires/makeewww 


^ IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in As¬ 
tronomy (AURA) under cooperative agreement with the National Science 
Foundation 
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only small changes in the abundances of the other elements. Bar¬ 
ium (scandium) were corrected by 0.05 (0.06) dex in UMI718, 
0.01 (0.03) dex in UMI396 and 0.02 (0.02) dex in UMI446. The 
copper abundance differed only in UMI718 (by 0.03 dex) with the 
inclusion of hyperfine splitting. 

Ap art from a Non-LTE cor rection applied to the Al/Fe abundance 
by ICohen & Huand (l2010h . all of the data used in this paper are 
based on the assumption of LTE conditions. Our stars do not show 
any signs of ionisation imbalance in the abundances we calcu- 
late for Fe and Ti w it hin th e uncertainti es. The calculations of 
iBergemann & Gehr^ (l2008h for Cr and iBergemann & Cescuttil 
(l2010h for Mn imply upward-corrections by ~0.4 dex in these el¬ 
ements, but we note that their parameter range mainly included 
dwarf stars and did not extend to cool giants as in our sample. Thus 
we refrain from applying these corrections to these elements in the 
following. 

The initial estimates for the temperature and the surface grav¬ 
ity listed in Table [T] ar e the photometric values ba sed on V-I, 
using the calibration of iRamirez & Mel endej (2QQ5h, wher e we 
adopted a reddening of E(B-V)=0.02 dSchlegel et al'l Il998h and 
an initial guess for the metallicity of [Fe/H] = -2.0. Through¬ 
out our analysis, we constructed model atmospheres for a large 
range of metallicities and microturbulent velocities by interpola¬ 
tion of Kurucz’^ grid of one-dimensional 72-layer, plane-parallel, 
line-blanketed models without convective overshoot, assuming lo¬ 
cal thermodynamic equilibrium (LTE) for all species. For the UMi 
stars, we inco rporated the opacity distr ibution functions ODFNEW 
described by dCastelli & Kurucz ll2003H^ . As UMI446 was found to 
be very metal poor, we tested our final results with the a enhanced 
AODFNEW models as well but found that these did not affect the 
results more than 0.01 dex for any of the elements. 

These model atmospheres were then used in MOOG, together with 
the EWs to calculate the abundances of the Fe I and Fe II lines for 
each star. The iteration to obtain the best model parameters was 
made first by checking the abundance trends with excitation poten¬ 
tial in order to estimate the spectroscopic temperature and remov¬ 
ing any trends with the reduced EW to estimate the microturbulence 
velocity. Finally, the best value for the surface gravity was mea¬ 
sured by minimising the [Fe I/Fe II] ratio. During this procedure, 
all steps were iterated to convergence. Furthermore, we excluded 
strong lines with EW>160 to avoid saturated lines, and EW<40mA 
to avoid biasses due to low S/N (e.g., Magain 1984), from this part 
of the analysis. Although it is difficult to calculate formal errors for 
the temperature, log g and microturbulence as these parameters are 
highly interdependent (though see Section [Jj]), our final models 
minimised the correlation coefficients for all of these (less than 6% 
for T and V for all of our stars), and we could achieve a reasonable 
ionization balance solution, where the difference between the abun¬ 
dances of Fe I and Fe II are 0.07, 0.05, and 0.03 dex for UMI718, 
UMI396 and UMI446 respectively. 

In a few cases where the number of measured lines were very 
small for an element, we included lines with EW<40mA if the 
lines were clear and without blends. These exceptional, weaker 
lines ([Ce/Fe] for UMI718, [Si/Fe], [Cr II/Fe II] and [Ce/Fe] 
for UMI396, [Cu/Fe], [Y/Fe] and [Nd/Fe] for UMI446) were re¬ 
checked during the iterative process of finding the stellar atmo¬ 
sphere parameters which were normally adjusted through the trends 
seen for the [Fe/H] abundances. Despite the very low S/N in 


Model 

Description 

Infall 

Galactic winds 

Final stellar 
mass [Mq] 

Mo [Mo] 

A 

Closed box 

No 

No 

~ 2.3 X 10^ 

1.2 X 10'' 

B 

Accretion 

Yes 

No 

~ 2 X 10^ 

2.4 X 10’ 

C 

Acc + winds 

Yes 

Yes 

~ 2 X 10^ 

2.4 X 10’ 


Table 4. The parameters in the chemical evolution models. Mo is the gas 
mass used in eq. [2| and the laws for gas infall and galactic winds are given 
in eq.[2|and eq.[3l respectively. 


UMI446, there were features strong enough to be measured in both 
blue and red parts of the spectra. 


2.3 Chemical abundance measurements and uncertainties 

By varying the model atmosphere parameters by their typical un¬ 
certainties of lOOK (T) and 0.2km/s (Vmic), we find the systematic 
errors in [Fe/H] to be 0.1 dex for UMI718, 0.15 dex for UMI718 
and 0.13 dex for UMI446. 

The abundances calculated for 17 elements in our observed UMi 
stars are presented in Table [3 together with t he So lar abundance 
scale, which was adopted from lAsplund et al.l (l2009h . The detailed 
comparison with previous data and our chemical evolution models 
are presented in Section [3] 

The overall abundance errors for each element given in Table [3 are 
determined by using cr/ ^^N where cr is the standard deviation of 
the individual line abundances and N is the number of observed 
lines for this el ement. In Table [3 we li st each of these values. 
Similarly to the lKirhv & roheJl2012h data, UMI446 and UMI718 
have radial velocities which make the telluric absorption lines co¬ 
incide exactly with the forbidden oxygen line at 630.03 nm. There¬ 
fore, this line was only measured in UMI396 by fitting synthetic 
spectra. 


3 ABUNDANCES AND THE CHEMICAL EVOLUTION 
MODELS FOR URSA MINOR 

In this section, we present the chemical abundances calculated for 
17 elements in three new UMi stars (see Table [3 and use a set of 
homogeneous chemical evolution models to investigate the impli¬ 
cations of the abundances in a data set comprising both our new 
data and those obtained in previous studies of UMi. The main char¬ 
acteristics of our models are summarised in Table|4] The difference 
among the models is the treatment of the gas flows. Model A is 
a closed box, without infall or outflow; model B has infall but no 
galactic winds; in model C, we take into account both infall and 
galactic win ds. We constra i n our models to follow the observational 
SF given bv ICarrera et al.l (l2002h . while imposing a star formation 
law given by 




[3 • 

-10“‘'’f + 0.5] 


0 


ift<5Gyr 
if t > 5Gyr 


( 1 ) 


where, is a normalized Gaussian function with cr of 0.5 Gyr, 
centered at to equal to 0.5 Gyr. The SF law provided by Eq.[T]niim- 
ics the observed peak of early SF with a subsequent period of ex¬ 
tended S F with lower efficie ncy, which stops after 5 Gyr as pro- 
www. stsci. edu/hst/observatory/crds/castelli_kurucz_atlas.htmlposed bv ICarrera et al.l (l2002h . 


http://kurucz.harvard.edu 
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SUN 


UMI718 



UMI396 



UMI446 


X 

10g6(X)o 

log6(X) 

cr(nr of lines) 

[X/Fe] 

10g6(X) 

cr(nr of lines) 

[X/Fe] 

log6(X) 

cr(nr of lines) 

[X/Fe] 

Fel 

7.5 

5.38 

0.17(88) 

-2.12+0.02 

5.37 

0.19(118) 

-2.13+0.02 

4.83 

0.26(68) 

-2.67+0.03 

Fe II 

7.5 

5.45 

0.19(5) 

-2.05+0.08 

5.42 

0.17(6) 

-2.08+0.07 

4.86 

0.19(7) 

-2.64+0.07 

O I synt 

8.69 

- 

- 

- 

6.99 

...(1) 

0.38 + 0.13 

- 

- 

- 

Mg I 

7.6 

5.74 

0.09(3) 

0.26 + 0.05 

5.79 

0.15(3) 

0.32 + 0.08 

5.26 

0.10(2) 

0.33 + 0.08 

Si I synt 

7.51 

- 

- 

- 

5.57 

...(1) 

0.19 + 0.13 

- 

- 

- 

Cal 

6.34 

4.55 

0.34(15) 

0.33 + 0.09 

4.35 

0.15(16) 

0.14 + 0.04 

3.97 

0.20(7) 

0.3 + 0.08 

Sc II* 

3.15 

1.14 

0.09(3) 

0.02 + 0.09 

1.08 

0.19(4) 

-0.02 + 0.12 

1.15 

0.27(3) 

0.63 + 0.18 

Til 

4.95 

2.96 

0.16(15) 

0.13 + 0.04 

2.86 

0.19(25) 

0.04 + 0.04 

2.73 

0.23(7) 

0.45 + 0.09 

Till 

4.95 

3.07 

0.09(15) 

0.17 + 0.08 

3.06 

0.15(17) 

0.19 + 0.08 

2.79 

0.34(13) 

0.48 + 0.12 

CrI 

5.64 

3.50 

0.18(12) 

-0.02 + 0.05 

3.27 

0.09(13) 

-0.24 + 0.03 

2.99 

0.07(6) 

0.02 + 0.04 

Cr II* 

5.64 

- 

- 

- 

3.33 

...(1) 

-0.23 + 0.13 

- 

- 

- 

Mn I* 

5.43 

3.16 

0.17(5) 

-0.15 + 0.08 

2.96 

0.12(3) 

-0.35 + 0.08 

3.05 

0.19(3) 

0.29 + 0.1 

Nil 

6.22 

4.24 

0.17(17) 

0.14 + 0.04 

4.11 

0.17(24) 

0.02 + 0.03 

3.87 

0.11(5) 

0.32 + 0.05 

Cu I** 

4.19 

1.92 

...(1) 

-0.18 + 0.1 

- 

- 

- 

1.64 

...(1) 

0.12 + 0.34 

Znl 

4.56 

2.41 

0.12(2) 

-0.03 + 0.08 

2.54 

...(1) 

0.11+0.13 

- 

- 

- 

YII 

2.21 

-0.32 

0.19(3) 

-0.48 + 0.13 

-0.35 

0.04(3) 

-0.48 ± 0.07 

-0.23 

...(1) 

0.18 + 0.34 

Ba II* 

2.19 

-0.46 

0.1(3) 

-0.55 + 0.1 

-0.35 

0.08(2) 

-0.45 ± 0.09 

-1.21 

0.47(2) 

-0.74 + 0.34 

Ce II* 

1.58 

-0.34 

0.03(2) 

0.13 + 0.1 

-0.65 

0.13(2) 

-0.15 + 0.12 

- 

- 

- 

Nd II 

1.42 

- 

- 

- 

-0.68 

0.14(3) 

-0.02 + 0.1 

-0.31 

...(1) 

0.89 + 0.34 

Sm II** 

0.96 

- 

- 

- 

-0.27 

0.04(3) 

0.9 + 0.07 

- 

- 

- 


Table 3. Abundances of chemical elements in our target UMi stars and the Solar abundance scale, which was adopted from lAsnlund et alJ i2009li .(*i For the 
values that come from the measurements of a single spectral line, we assume the abundance error to be the same as the largest calculated error for that star 
(0.1 for UMI718 and UMI396, 0.3 for UMI446). The abundance values are given in columns 5, 8 and 11 in terms of [X/Fe] and [Fe/H] for the Fe I and 
Fe II. (**)We measure a total Sm II abundance of [Sm/Fe]=0.88 for UMI396, with a surprisingly low standard deviation of 0.04 dex from the lines at 451.96, 
452.39 and 456.62 nm. Although interesting, as all three lines are very weak with EW<25 and are at the blue end of the spectra, we omit this element in our 
interpretation of the chemical abundances in UMi. 


The infall law used in models B and C takes the same form as the 
peak of the SF described in Eq[T]and is given by 

Ginfaii(0[Mo/yr] = M, • (2) 


Models B and C have no initial gas and reach Mo which is given 
in Table |4] following the infall law given in Eq. |2] Although the 
infall law is not directly included in the SE law, the common form 
of these equations provides an implicit connection between the SE 
and the inf all. 

The galactic wind for model C has the simple form 


Gwi„d(f)[10-Wo/yr] 


0 

2 • 10-5 . 


ift<2Gyr 

ift>2Gyr 


and is assumed to act only at the end of the gas infall and of the 
main SE episode (after 2 Gyrs). By implementing the winds with 
the simplest possible assumptions, we are able to investigate their 
implication on the chemical evolution, without being forced to con¬ 
strain their origin(s), given the fact their origins are still not clearly 
understood. 

The top row in Eig.[T]shows the total mass, gas mass and the stellar 
mass for the three models as a function of time, while in the second 
row we show the three SE histories obtained for each individual 
model. As model A, which has all the gas available from the outset, 
produces stars more efficiently than models B and C, in order to 
obtain the same final stellar masses in all models, the initial gas 
mass in models B and C is twice that of model A. Both models A 
and B have a large amount of gas left at the end, however, and it is 
clear that the loss of gas through winds is necessary to explain the 
lack of gas observed in UMi. 

In all of our models the IME of ISalpetei] (Il955h is used, and 
the other chemical evolution parameters, such as the treatment of 


SNe la, and the n ucleosynthesis, as well as the stellar lifetimes are 
those adopted by lErancois et al.l (l2004l) . The excep tions are Mn, 
for wh ich we follow the nucleosynthesis adopted in ICescutti et al.l 
(l2008h . and the neutron capture elements Ba and Y, for which 
we adopt the same nucleos ynthesis of model EC+s 2 shown in 
ICescutti & Chiappinil (l2014l) . The latter study assumed that the r- 
process elements are produced by electron capture SN in stars with 
masses between 8 - lOM© and the s-process elements are pro¬ 
duced in rotating massive stars (spinstars) with masses between 
15and40Mo. 


The parameters we assumed aim to reproduce the observed metal - 
licity distribution function (MDE) observed bv iKirbv et al.l (l201lh . 
Eig.|2] compares the MDE found in our models with the observed 
distribution of stellar metallicities . By m eans of a simplified chem¬ 
ical evolution model, iKirbv et ^ (l20Ilh already showed that it is 
not possible to find a model which is able to reproduce the observed 
MDE without infall. Our findings agree with their conclusions, and 
indeed our model A only poorly approximates the observed MDE. 
While our model B (which includes gas infall) is able to match the 
main characteristics of the observed MDE, as shown by Eig. [T] it 
cannot explain the absence of gas displayed at the present time by 
UMi - which is also a common feature of all the dSph galaxies 
around the Milky Way. The implementation of a galactic wind in 
Model C, albeit in a simplified manner, solves this problem. We 
note that it also improves the match of the theoretical MDE with 
the observed one by slightly extending the tail of the distribution 
at high metallicity. We emphasize that the MDE and the SE law 
are the only observational constraints that were used in setting the 
constants in Eqs. [Tl2l and 0 It is then noteworthy that the ensue- 
ing chemical evolution model, shown in Eigs. [3ll7l are successful in 
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Figure 1. Top: The mass of baryons, gas, stars and wind as functions of time for models A, B and C (from left to right). Bottom: the SF rate as a function of 
time for the same models. 




S e 2. Comparison between the observed MDF histogram (iKirbv et al.l 
and the theoretical MDF produced by our models A (black), B(blue) 
and C(red). 


Figure 4. Oxygen abundances in 12 UMi stars . Literature data are 
taken from IShetro ne et al.l (l200lh (black arrows for the upper limits) and 
ICohen & Huand (l201C)h (green dots). Our star is shown with the red hlled 
circle. The chemical evolution models are the closed box (black dotted line), 
model with 5 Gyr of gas infall (blue solid line) and the model with infall 
and gas loss through winds (red dashed line). 


matching the observed abundance ratios, even though they were not 
used in fitting the model parameters. 


3.1 Qf-elements: O, Mg, Si, Ti, Ca 

The measured cr-element abundances (see Fig.[3]) are generally con¬ 
sistent with a SF efficiency that is lower than that of the Galac¬ 
tic halo, as also shown by our models which show a gradual de- 
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[Fe/H] [Fe/H] 


Fig ure 3. a-elements in UM i. The symbo l are data from : this w ork (red hlled circles). Ishetrone et ~S1 i200l| ) (black squ ares, and arrows for the upper limits) 
and ICohen & Huanj i201(j) (green dots). ISadakane et aP 1 2004|) (empty mag enta diamonds) and iKirbv & Cohei] i2Q12h (blue empty hexagon). The average 
abundances for these elements in halo stars as calculated by Cavrel et ai] <2004 are: [Mg/Fe]=0.26 dex, [Si/Fe]= 0.42 dex, [Ca/Fe =0.35 dex, [Ti/Fe]=0.26 dex. 


crease in the [cr/Fe] caused by the onset of SNIa already at very 
low metallicities (in the models at [Fe/H]= -3). On the other 
hand, the data for both [Mg/Fe] and [Ca/Fe] indicate a plateau at 
[Fe/H]< -2 dex though sparsely sampled, which is not found in 
the models. For the rest of the individual elements, our [Ca/Fe] 
and [Ti/Fe] abundances, w hich are closer to lShetrone et al.l (I 2 OOII) 
and ISadakane et aL (l2004h data at higher [X/Fe] values than the 
ICohen & Huand feOld) data, are only slightly lower than the halo 
abundances. The [Si/Fe] ratio has a larger scatter in UMi than in 
the halo (where it is nearly constant), and the [0/Fe] ratio which is 
reproduced well by our models (see Fig, is lower in UM i than 
in the halo stars dCavrel et al.ll2004l : ICohen & HuandbOld) . How¬ 
ever, additional lower metallicity stars are needed in order to better 
distinguish among the different models. 


As seen in Fig.[3l model A (without infall and winds), tends to have 
a less pronounced decrease at higher [Fe/H] which is due to the 
relatively higher SF compared to models B and C. Concerning these 
two models, the only difference is that model C (with winds) is able 
to produce stars at a higher metallicity. This reflects the MDF we 
obtain for it, driven by the late metal pollution (mostly SN la) in a 
decreasing reservoir of gas which produces a higher concentration 


of metals. On the other hand, we underline that while the winds do 
not change significantly the chemical evolution for this galaxy, they 
help explain the absence of gas at the present time and improve the 
match with the MDF. Instead of evoking an episodic star formation, 
our models use an extended, low-efficiency star formation period. 


3.2 Fe-peak elements: Sc, Cr, Ni, Mn, Cu, Zn 

As seen in Fig.[5l our models generally reproduce the abundances 
of the Fe-peak elements well. Similarly to the ISadakane et al.l 
(I 2 OO 4 I) data, our siSLVs hav e abundances close to those found in halo 
giants dCavrel et al.ll2004h in Ni and Cr (slightly highe r in Zn ), and 
are higher compared to those from ICohen & Huan3 (I 2 OIC 1 I) . The 
[Cr/Fe] abun dance at nearly the so lar value in UMI446, UMI718, 
UM I K from Shetrone et al.l (I 2 OOII) as well as the two more metal 
rich ISadakane et al.l (l2004l) stars shows a plateau for this element 
apart from the decline towards lower metall icities which is likely to 
be cau sed by NLTE conditions as argued bv iBergemann & Cescuttil 
(I 2 OIQI) . The Cu abundances we measure in UMI718 and UMI396 
are in agreement with the whole UMi sample available to date, 
which show a decline with increasing metallicity. 































































Chemical abundances and the evolution of Ursa Minor 1 


1.0 


0.5 


St 0.0 
tr 


-0.5 


- 1.0 



k 

. I .s 


XX .. ‘ , X 

. closed box 


— accretion 


■ -- acc+winds 

_1_1_1_1_ 


-4.0 -3.5 


-3.0 -2.5 

[Fe/H] 


-2.0 -1.5 -1.0 


1.0 


0.5 



- 1.0 

-4.0 -3.5 -3.0 -2.5 -2.0 -1.5 -1.0 

[Fe/H] 






Figure 5. Abundances for the Fe-peak elements. The symbols and colors are the same as in Fig. [3] 


For Mn, we match the bulk of the data at intermediate metallic- 
ity, so the model is relatively successful, however the decrease pre¬ 
dicted at the extremely low metallicities by our models is not con¬ 
sistent with the data in which this signature is detected at [Fe/H]=- 
3 dex. Our most metal poor star, UMI446, deviates from rest of 


the data set with a strong enhancement in Fe-peak e bments, partic¬ 
ularly in Sc and Mn. Although stars measured by ISadakane et al.l 
(I 2 OO 4 I) display [Mn/Fe]>0 as well, UMI446 has a much lower 
[Fe/H] than these two, and cannot be explained by our chemical 
evolution model even if we take into account metal-dependent Mn 
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Figure 7. The ratio of yttrium to barium abundances in the UMi stars. The 
symbols and the colors are as in Fig. [3] 


yiel ds for both SNe la and SN II as used by ICescutti et al.l (l2008h 
an d lNorthetal.lJ2012h . 

As a final remark in this subsection, we find the chemistry shown 
by the most metal rich star measured by ICohen & Huand (1201 Oh 
challenging. An anomalous pollution by SN la which cannot be 
simulated by our simple chemical evolution models, could explain 
the measured drop in the ratio of cr-elements (and neutron capture 
elements). However, it is more difficult to explain the strong enrich¬ 
ment of Fe in this star without an enrichment of the other Fe-peak 
elements. 


3.3 Neutron capture elements: Y, Ba, Ce, Nd 

There are several interesting points that are seen in the neutron cap¬ 
ture element abundances in UMi. First, all of the neutron capture el¬ 
ements that we measure have very large scatter across the observed 
stars, covering a range of about 2 dex. The Y and Ba abundances 
for the data set are shown in Fig. [6] and the [Ce/Fe] and [Nd/Fe] 
abundances in two new stars each are given in Table [3] Due to this 
spread in the observed abundances, we expect that inhomogeneous 
modeling would be required to obtain a more precise match to the 
data (see for example Cescutti 2008, Cescutti et al. 2013 and Ces- 
cutti and Chiappini 2014). Nevertheless, a simple chemical evolu¬ 
tion model should be able to reproduce the mean trend. 

The second point of note is the increase of [Ba/Fe] with metallicity 
seen in the entire UMi sample (Fig.[6l) and which o verlaps with the 
halo stars (lFulbrightll2000l : l^hen & Huand l2010h . In Fig. [6l we 
show that the model results for Y and Ba are in reasonable agree¬ 
ment with the data in that there is a systematically lower average 
[Ba/Fe] ratio at [Fe/H]«3 dex compared to higher metallicities. In 
addition, although the present data in UMi cannot resolve these de¬ 
tails, the models predict a dip in both the [Y/Fe] and [Ba/Fe] abun¬ 
dance ratios at [Fe/H]-3, which is due to the time scale at which 

the r-p rocess production by electron capture SN (^see lCescutti et al.l 
(1201^ for details) starts. 

Finally, as predicted by the theoretical nucleosynthesis prescrip¬ 
tions we use, the data show a much lower [Ba/Fe] ratio at very low 
metallicity, where the [Y/Fe] remains approximately at the same 
level. At lower metallicity, the production of these neutron capture 
elements could be sustained by s-process production in fast rotat¬ 


ing stars. Note that the low [Ba/Fe] abundance ra tio is also found in 
the Hercules dSph over a broad range in [Fe/H] (iKoch et al.ll2013h 
This effect is clearer when we plot the ratio between these two 
neutron capture elements (see Fig.[71). In this model, before the r- 

process production starts at [Fe/H]-3, the ISM is enriched only 

through s-process production in spinstars, which is confirmed by 
an overabundance in [Y/Ba] found in the data at very low metallic¬ 
ities. There is possibly an offset of ~0.5 dex between the last data 
with high [Y/Ba] defined b y UMI446 from this work and COS4 
from ISadakane et aP (l2004l) and the knee found in the model; this 
can be explained both in terms of inhomogeneity of the system and 
also in tiny variations in the parameters of the chemical evolution 
model, definitely inside the uncertainties of our modelling. 


4 DISCUSSION 


In this paper, we presented a new data set of chemical abundances 
for the Ursa Minor (UMi) dSph, based on spectra taken with the 
Keck HIRES spectrograph, as well as a new set of chemical evolu¬ 
tion models for 12 elements. A data set combining our new mea¬ 
surements with data from previous studies is consistent with a sce¬ 
nario in which the very metal poor stars in the dSphs have a- 
element abundances generally comparable to those in the Galactic 
halo, although in UMi, there is a gradual decrease in the a elements 
with increasing metallicity which indicates that the ISM was also 
enriched by SN la. The halo stars do not show this trend which is 
a signature that the star formation in the halo stopped on a shorter 
timescale compared to the timescale for the enrichment by SN la. 
Our chemical evoluti on models use s i milar nucleoynthesis pre¬ 
scriptions to those by lUanfranchi et al.l (l2006l. l2008h for modeling 
UMi with the smaller data set available at the time. However, our 
implementation of the star formation rate is dif ferent as we derive 
the SFH exclusively from the CMD, as done bv iNorth et al.l (l2012h 
for computing the manganese abundances in four other dSphs. 
Here, we extended these latter models to study Ursa Minor for the 
first time and to include the other elements for which we have mea¬ 
sured abundances. The most striking novelty in the models is the 
nucleosynthesis adopted for the neutron capture elements; due to 
the lack of data, the signature of increasing [Y/Ba] ratio toward 
lower metallicities was not studied in these previous works. The 
overabundance of [Y/Ba] at very low metallicity, similar to that 
found in the Galactic halo, is reproduced by nucle osynthesis in the 
spinstars as argued in ICescutti & Chiappinil (l2014h and is modeled 
for the first time for a dSph. This trend is due to the timescale of 
the onset of r-process prod uction by electron capture in SN (see 
ICescutti & ChiaDDinill2014l) : before this point the production of the 
neutron capture elements could be sustained by s-process produc¬ 
tion in fast rotating stars dCescutti et al.ll20I3b . 

In the ratio of the [cr/Fe] abundances, it is not easy to identify the 
trends in the data, ie. whether these elements show plateaux at low 
metallicities up to around [Fe/H]= -2 dex, or whether they are 
matched better by a gradual decrease of [cr/Fe] starting from very 
low metallicities as found in our models. In either scenario, together 
with subsolar values of [Ba/Fe] up to the same metallicity where it 
reaches the solar value, the abundances are consistent with a SF 
process that lasted a relatively long time (^ 5 Gyrs) but had a very 
low efficiency until it stopped at [Fe/H]=-1.3 dex. While earlier 
studies found that an inefficient SF in a system without infall or 
mass loss seemed to fit the data dlkuta & Arimoto 2002h . as the 
spectroscopic samples grew, already [Kirby et al. d2013h were able 
to show that to obtain the MDF, a model with the infall of gas was 
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Figure 6. Abundances for the neutron capture elements in UMi giants.The symbols and colors are the same as in Fig. [3] 


required. Our best model is one where a gas loss proportional to the 
amount of total gas was also included (winds), and could account 
for both the low efficiency of the SF at later stages and the observed 
lack of gas today. We emphasize that, realistically, the gas loss from 
the dSph is due to a combination of supernova feedback, winds and 
tidal interactions and hence would have a more drastic effect. 

Our results are consistent with the previous data except for a devia¬ 
tion of the abundances of the Fe-peak elements for UMI446 which 
is our most metal poor star. The [Fe/H] of this star is very simi¬ 
lar to COS4 from ISadakane et all (l2004f) which has [Fe/H]= -2.7 
dex. Although abundances of the low Z cr-elements are consistent 
with the rest of the sample, UMI446 has an unexpectedly higher 
abundance of Sc and the of-elements Ti and the Fe-peak elements 
Cr, Mn, Ni, Cu, Zn. At Z=30 it becomes consistent again with the 
rest of the data, which is especially important in the relation of the 
[Y/Ba] ratio with metallicity. Among these, one of the most chal¬ 
lenging features is the overabundance of Mn for UMI446. These in¬ 
triguing results could point to an inhomogeneous chemical enrich¬ 
ment of the ISM at the early stages. Although this inhomogeneity 
would be short lived, it might be able to produce a star like UMI446 
if the explosive nucleosynthesis of a SN II through a very massive 
progenitor wi th bipolar explosion (iet-li ke) is taken into account, as 
suggested by iMaeda & Nomotol (l2003h . 

In the future, we plan to investigate the production of the spread in 
the observed abundances in UMi using more sophisticated, inho¬ 
mogeneous models to have a clearer insight to its chemical evolu¬ 
tion. 
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